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Various forms of matter may be produced in ultra-relativistic heavy ion collisions. These are the Quark Gluon 
Plasma, the Color Glass Condensate , the Glasma and Quarkyonic Matter. A novel effect that may be associated 
with topological charge fluctuations is the Chiral Magnetic Effect. I explain these concepts and explain how 
they may be seen in ultra-relativistic heavy ion collisions. 



1. Introduction 

The central purpose of ultra-relativistic nuclear col- 
lisions is to produce and study forms of high energy 
density strongly interacting matter and to study their 
properties. A variety of such forms of matter have 
been proposed. (References to early developments in 
this field may be found in Ref. [T]) The Quark Gluon 
Plasma (QGP) is a deconfined phase of quarks and 
gluons [2 . The Color Glass Condensate (CGC) is high 
energy density and highly coherent matter made from 
gluons [3]. The Glasma is matter produced in the 
first instants of a heavy ion collision from the Color 
Glass Condensate [l]-[9]. It involves many strongly in- 
teracting and highly coherent color electric and color 
magnetic longitudinal flux lines. Quarkyonic matter 
is high baryon density matter where the energy den- 
sity is very high compared to the typical QCD scale, 
but that still remains confined }93|-|94], I will describe 
each of these forms of matter and how they might 
appear in ultra-relativistic nuclear collisions. 

A simple picture of ultra-relativistic nuclear colli- 
sions is shown in Fig. [T]- Here two Lorentz contracted 
nuclei, which arc thought of as sheets of Color Glass 
collide. In the instant these two sheets pass through 
one another, the fields become singular, and produce 
the Glasma. The Glasma fields decay and produce a 
QGP. This QGP dilutes through expansion and even- 
tually form a gas composed of ordinary strongly in- 
teracting mesons and nuclcons. 



tions are associated with the generation of masses for 
the nucleons. Analogous processes in electroweak the- 
ory may be responsible for the baryon asymmetry of 
the universe. 

There are many theoretical questions that we would 
like to answer through the collisions of such matter. 
Among them are: 

• What is the high energy limit of QCD? 

• What are the possible forms of high energy den- 
sity matter? 

• How do quarks and gluons originate is strongly 
interacting particles? 
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Figure 2: A lightcone diagram for ultra-relativistic 
hadronic collisions. 
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Figure 1: Ultra-relativsitc nuclear collisions and various 
forms of high energy density matter. 

Sometimes it is useful to conceptualize such colli- 
sions in terms of a light cone diagram as shown in Fig. 
[2] This figure shows the various times in the collision. 
In particular, in the Glasma phase, color electric and 
color magnetic fields produced can carry topological 
charge. [T2]-[TT] In QCD, topological charge fluctua- 



2. The Color Glass Condensate 

A remarkable property of strongly interacting par- 
ticles is that the small x part of their wavefunction is 
dominated by gluons. Since x is the ratio of the energy 
of a constituent to that of the hadron in the reference 
frame where the hadron is fast moving, the smallest 
values of x are probed at the highest energies. It is 
this gluon rich part of the hadron wavefunction that 
controls the high energy limit of QCD. The distribu- 
tion of gluons as a function of x at Q 2 = 10 GeV 2 is 
shown in Fig. [3] for a proton. 

Because the size of a proton grows at most like 
ln 2 (E/E ) at high energies, and because the gluonic 
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Figure 3: The gluon and quark composition of a proton. 



This high energy density saturated gluonic matter 
is called the Color Glass Condensate t 20]- [26]. The 
word color arises from the color of the gluons. The 
word glass is because the coherent fields at small x 
are ultimately generated by gluons at larger values of 
x through evolution equations. These fast gluons have 
their natural time scale dilated, and this dilated time 
scale is then passed on to the fields at small x. Sys- 
tems of particles that have long time scales for evolu- 
tion compared to natural time scales are glasses. The 
word condensate is because the phase space density of 
gluons is very large and they gluons are very coherent. 
An effective energy functional in terms of the phase 
space density, 



dN 



dyd 2 pTd 2 XT 



(1) 



contribution to a hadron is measured to grow much 
faster, like a power l/x & where 6 = 0.2 - 0.3, the glu- 
ons must form a high density state. The coupling 
evaluated at such a high density scale is therefore 
weak. Weak coupling does not however mean weak 
interactions. Coherence can amplify the effects of an 
intrinsically weak interaction, as we know is the case 
in gravity. 

Surely if we try to insert more and more gluons of 
fixed size into a hadron wavefunction, at some point 
it becomes very difficult to add more |T5]-[TS] . This 
is when the interactions of order as become strong 
because there are typically 1 /as particles at that size 
scale. This requires that additional particles added 
to the system are at yet a smaller size scale. This is 
shown in Fig. [4] We see that saturation requires that 
the number of gluons of fixed size will stop rapidly 
growing above some energy, but the total number of 
saturated gluons grows forever, since one saturates 
gluons of smaller an smaller size. 



is 
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Figure 4: Saturation of gluons in the hadron wavefunction. 



E 



-Kp + k asp 



(2) 



The first term in this equation drives the condensa- 
tion and the second is due to repulsive interactions 
that stabilize it. The extrcmum is when p ~ l/ccg. 
The dependence on coupling is typical of condensa- 
tion phenomenon such as Bose condensation, super- 
conductivity, and even the Higgs condensate. 

It may seem a little strange to think of a hadron 
at high energy as a sheet of Color Glass Condensate 
and at low energies as valence quarks. It naively vio- 
lates Lorentz invariance. The resolution of this appar- 
ent paradox is that these correspond to two separate 
pieces of the hadron wavefunction. The matrix ele- 
ments important for high energy scattering are sen- 
sitive to the high gluon density part of the hadron 
wavefunction, where low energy processes are sensi- 
tive to the valence quark piece. 

The typical density of gluons that arise from the 
phase space density above is 



dN 
dy 



= / d 2 x T 
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where Q sa t is the maximal momenta for the gluons 
that are fully saturated, that is where p ~ l/as- 



3. The Initial Conditions and the Glasma 

The color electric and color magnetic fields in 
a sheet of Colored Glass are the Lorentz boosted 
Coulombic fields, Their color and intensity is random, 
and they are plane polarized perpendicular to the di- 
rection of motion of the hadron. They are analogous 
to the Lienard-Weichart potentials of electrodynam- 
ics. The collisions of two nuclei may be thought of as 
the collision of two sheets of colored glass as shown 
in Fig. [5] Note that color electric and color magnetic 
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Because there are both longitudinal color electric 
and longitudinal color magnetic fields, there is a non- 
zero density of topological charge, E ■ B. A non-zero 
average of such a quantity would signal the breaking 
of P and CP. Here the average over space is zero, al- 
though there are strong local fluctuations. A topolog- 
ical charge generates non-zero helicity in quarks. To 
see this, note that if we accelerate a charge along the 
direction of E, it will spiral with a definite handed- 
ness around B. If we reverse the charge, the direction 
reverses and the spiral reverses so the helicity remains 
the same. 



4. The Quark Gluon Plasma 



Figure 5: The collision of two sheets of colored glass. 



fields are treated on an equal footing, so that the CGC 
is self-dual under E <-> B 

If one takes the initial conditions for the CGC fields, 
and then computes the field infinitesimally in time af- 
ter the sheets of color glass pas through one another, 
one finds longitudinal color electric and color mag- 
netic fields between the two sheets. It is as if the two 
sheets of colored glass became dusted with equal and 
opposite densities of color electric and color magnetic 
charges. We can think of the lines of color electric and 
color magnetic flux that join these sheets as flux tubes. 
The typical transverse size of these random flux tubes 
is R ~ l/Qsat- Such flux tubes are shown in Fig. [6] 
This collection of flux tubes and its subsequent evolu- 
tion is called the Glasma, since it is a transient state of 
matter that connects the Color Glass Condensate to 
the Quark Gluon Plasma. It has coherent fields, like 
the CGC, but it has properties similar to a plasma 
with strong background fields, and these fields evlove 
in time with the natural time scale 1/Q sa t- 




Through some as yet not completely understood 
mechanism, the Glasma eventually thermalizes into 
a Quark Gluon Plasma. [27] The QGP is matter that 
is deconfined and chirally symmetric 28J . There is ev- 
idence that this matter is well thermalized, since the 
distributions of produced particles are well described 
by hydrodynamic equations with surprisingly small 
viscosity [2], [25], [3D]- [33]. This is only understood if 
the constituents of the QGP are strongly interacting. 
A plot that shows the typical times scales and energy 
densities believed to be present in RHIC collisions of 
gold nuclei is shown in Fig. [7] Note that the maximal 
energy density is many times that of the typical scale 
of QCD e~ 1 GeV/Fm 3 . 
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Figure 7: Time scales for heavy ion collisions. 



5. Phenomenological Evidence for the 
Color Glass Condensate 



Figure 6: Glasma flux tubes. 



The theory of the Color Glass Condensate provides 
a framework for understanding the high energy limit 
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of QCD. This framework can be proven from first prin- 
ciples in QCD if the saturation is sufficiently large so 
that weak coupling methods can be used. A substan- 
tial question is: When is this true? 

There are a number of pieces of evidence that 
come from the Hera experiments on deep inelas- 
tic scattering, and from fixed target experiments on 
nuclear targets.. |34j- [35] Describing these results in- 
volves knowing the deep inelastic structure function 
F%, structure functions for diffractive deep inelastic 
scattering, and the longitudinal structure function Fl ■ 
A simple comprehensive description of these data are 
possible within the CGC framework [5rJj . 

One measure of the success of the CGC descrip- 
tion is its simple prediction of geometric scaling, that 
is that the cross section for virtual photon scatter- 
ing from a protons should scale as F(Q 2 /Q 2 aat )(x) 
|37j-|41j. The saturation momentum's dependence 
on x has been computed and agrees with the fits to 
the Hera data [H]-j33]. A geometric scaling curve is 
shown in Fig. [8] 



laboration, shown in Fig. [9] [45 . 
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Figure 8: Geometric Scaling in Deep Inelastic Scattering 

In spite of the remarkably good and simple descrip- 
tion that the theory of the CGC gives for the Hera 
data, there is no consensus within this community 
that this description is the correct one. 

One of the central results of the theoretical descrip- 
tion of the CGC is the existence of a renormalization 
group description as a function of x or rapidity. |44j 
An implication is that there should be an approximate 
limiting fragmentation or extended scaling. Such ex- 
tended scaling has been observed by the Phobos col- 



Figure 9: Extended scaling in deep inelastic scattering 

One of the earliest successes of the CGC description 
of heavy ion collisions was its predictions for the total 
multiplicity of produced particles and its dependence 
upon the centrality of the collision [6], [46] -[47]. Be- 
cause the saturation momentum provides an infrared 
cutoff in the computation of particle production, the 
total multiplicity is reliably computed if the satura- 
tion momentum is sufficiently large. It predicts a slow 
growth with energy and with centrality of the collision 
for the number of particle produced per nucleon par- 
ticipant in the collision. 

The CGC gives a heuristic derivation of the slow 
logarithmic growth of total hadronic cross sections 
[IB]- [39]. It provides a theory of non- leading twist 
shadowing, where the saturation momentum appears 
as the scale factor for the non-leading twist contribu- 
tions [50 - [54 . This theory was successfully applied to 
and predicted the suppression of particles produced at 
small nuclear x in deuteron gold scattering [55]- [56] . 
It will be severely tested when data is available from 
the recent dAu run at RHIC for its predictions con- 
cerning J /ty productions and forward backward an- 
gular correlations for particles in the small x region of 
the gold nucleus . The prediction for dAu J/ 1 ^ data 
is that when the saturation momentum exceeds the 
charm quark mass, then charm is like a light mass 
hadron. It should have approximate limiting frag- 
mentation, and production cross sections should be 
reduced relative to incoherent scattering. [57] For two 
particle correlations at small nuclear x , if one trig- 
gers on a high transverse momentum particle [58 -[59 , 
then there should be broadening and loss of intensity 
of the backward going peak relative to pp collisions 
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due to the interactions during the production of the 
backward going particle. Leading twist shadowing, 
that is perturbative QCD, predicts no such suppres- 
sion. 



6. Predictions of the Glasma 



emitters. A sum of negative binomial distributions 
is a negative binomial distribution. The Phenix ex- 
periments has measured the fluctuation spectrum of 
produced particle in heavy ion collisions and it is nega- 
tive binomial. [BS] The number of emitters scale as the 
number of participants, as expected for the Glasma. 



An important feature of the Glasma is long range 
correlations in rapidity induced by the Glasma flux 
tubes.|60]-|63). Correlation over a longitudinal mo- 
mentum scale p must be generated at a time t ~ l/p, 
a statement that can be proved using general causality 
arguments. In two particle correlation experiments at 
RHIC, a long range rapidity structure has been ob- 
served [55] • The plot is in pseudo-rapidity and az- 
imuthal angle. It is seen both in two particle correla- 
tions where one particle is triggered at high momen- 
tum, and inclusively, where there is no minimum mo- 



mentum requirement. In Fig. 10 preliminary results 
on the inclusive ridge from STAR is shown. 

The data show that for peripheral collisions the two 
particle correlation is similar to that in pp, where 
there is no long range component in rapidity. As the 
centrality increases a component extended over sev- 
eral units of rapidity appears. This component has a 
particle composition similar to single particle inclusive 
measurements. The STAR experiment has a limited 
rapidity coverage. Over a much longer range in rapid- 
ity, the Phobos experiment shows a ridge, as shown in 
Fig. " 
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The ridge can arise from an underlying flux tube 
structure. The long range correlation in rapidity is 
built into the longitudinal color electric and color mag- 
netic flux tubes. The azimuthal angular correlation 
can arise form the flux tube flowing with the un- 
derlying media in which it was produced. The flow 
Lorentz boosts the distribution and narrows an az- 
imuthally symmetric angular distribution in the di- 
rection of overall motion of the flux tube. 

The STAR experiment has measured the strength of 
long range correlations in rapidity. |66J These increase 
with the centrality of the collision and for peripheral 
collisions agree with pp collisions. For central gold 
gold collisions at the highest energy, the correlation is 
as strong as it can be based and is three times stronger 
than the upper limit predicted by impact parameter 
fluctuations. The dependence upon centrality can be 
understood as trade off of flux tube emission from the 
Glasma that is leading order in as and short range 
scattering. In more central collisions or higher energy, 
as is smaller, leading to a relative strengthening of 
flux tube emission and hence a stronger correlation 

EZI-EB1- 

The distribution of quanta produced from a single 
flux tube is a negative binomial distribution. A neg- 
ative binomial distribution is parameterized by the 
number of particles per emitter and the number of 



7. Topological Charge and the Chiral 
Magnetic Effect 

The fields in the Glasma carry a topological charge 
density since there are parallel color electric and color 
magnetic fields [H]-[T7]. Such a distribution has to- 
tal zero topological charge, but even total topological 
charge may be generated by evolution of the mater dis- 
tribution in the Glasma from fixed initial conditions. 
Conifigurations with topological charge may generate 
net helicity for quarks. In an overall net (electromag- 
netic) magnetic field, the magnetic moments of the 
quarks will align. If there is net helicity, then there is 
net motion of the quarks move along the direction of 
their spins, and an electromagnetic current is set up. 
This is the Chiral Magnetic Effect. 

In off impact parameter zero heavy ion collisions a 
net magnetic field is set up in the region where the nu- 
clei overlap. Such a field decays away quickly in time 
and can be computed from knowledge of the distribu- 
tions of ahcarged particle produced in such collisions. 
The result of a computations of this field is shown in 
Fig. |12| It is quite intense at early times and can have 
a characteristi scale of order 100 MeV. If there is a 
topological charge fluctuation that generates a net he- 
licity, this will therefore appear as a fluctuation in the 
net charge of particles perpendicular to the reaction 
plane. 

Of course any effect in QCD must be P and CP 
odd, so effects generated by the chiral magnetic effect 
must average to zero. Nevertheless, their fluctuation 
spectra nay be predicted by QCD, but they must be 
carefully compared against other possible effects. One 
measurable is the fluctuations in particle perpendicu- 
lar to the reaction plane. This has been measured in 



STAR and is shown in Fig. [13] [71]. 

The results shown agree semi-quantitatively with 
predictions of the chiral magnetic effect 



8. The Decay of the Glasma and 
Thermalization 

The Glasma configuration will decay away into glu- 
onic quanta. This decay is chaotic and turbulent |72j - 
|73j Two field configurations initially very similar de- 
viate from one another exponentially in time. This is 
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Figure 10: The ridge as seen in the STAR experiment 



Au+Au 200 GeV, - 30% 




Figure 11: The ridge as seen in the Phobos experiment. 



due to modes dependent on rapidity. The lowest order 
Glasma solution is taken to be invariant in rapidity, 
but the fluctuations will generate deviations from this. 
This has led to a number of speculations about turbu- 
lent thcrmalization Such thermalization might lead to 
anomalously low viscosities as suggested by the RHIC 
data. It might also lead to significant jet energy loss 
at early times as jets traverse the Glasma fT4"]-|76|. 

The problem with all of these speculations is that 
the Glasma lives a very short time t ~ 1/asQsat be- 
fore the effects of any coherent field should be drained 
away. On the other hand, thermalization involves 
scattering, and we would expect time scales of roder 
t ~ l/a|Qsat for this to occur. If the matter is ther- 
malized, it has little coherent component, and the de- 
scription of distributions of particles and flow using 




0.5 1 l.fl J 2.5 3 

T(flll) 



Figure 12: The typical magnetic field strength produced 
in the overlap region in the collision of two Gold nuclei. 

hydrodynamic simulation suggests that such matter 
has very small viscosity, and hence is very strongly 
interacting. 

Another possibility is that the thcrmalization is a 
consequence of the dynamics becoming very strongly 
coupled as the matter expands. This strongly coupled 
quark gluon plasma may have qualitative features in 
common with that of the AdSCFT model of QCD 
interactions. This theory is TV = 4 supersymmetric 
Yang-Mills theory and as such has a peripheral re- 
lation to QCD. The theory predicts a lower bound 
on the viscosity to entropy ratio that is a factor of 
2-3 times smaller than that which seems to describe 
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experimental data within hydrodynamic simulation. 
[77] The lower limit found in AfSCFT computations 
is close to that found by uncertainty principle bounds, 
that is, the mean free path must be larger than the 
Compton wavelength of a particle. 

AdSCFT computations have proved useful in a 
more theoretical context: The naive extrapolation of 
perfect fluid hydrodynamics to include viscous terms 
violates either causality, or the requirement that en- 
tropy is always increasing as a function of time. The 
AdSCFT computations suggested specific forms for 
higher order modifications of the hydrodynamic equa- 
tions that do not violate these principles, and provide 
estimates for the magnitude of coefficients that char- 
acterize these corrections [78]-[8"2]. 

One of the outstanding problems in the description 
of the QGP is a quantitative evaluation of jet energy 
loss. [83] Jets are suppressed by a factor of 4-5 in cen- 
tral Au-Au collisions at RHIC energies out to an en- 
ergy of at lest 15 GeV. Heavy quarks such as charm, 
are suppressed to the same degree as light quarks, 
and appear to flow with the media produced like light 
quarks. These phenomena are difficult but possible to 
include in weak coupling QCD computations. 

Perhaps a more natural explanation arises in the 
AdSCFT strong coupling world. The problem is that 
in this strongly couple world, jets are produced neither 
in electron-positron annihilation nor in deep inelastic 
scattering. If one tries to produce a jet, it immediately 
produces particle and cascades down to low momen- 
tum, producing isotropic distributions of low momen- 
tum quanta. This has the consequence that the dis- 
tribution of parton in a hadron scale as 1 /x 2 , and all 



the energy in a hadron is in the low x quanta. In fact 
these quanta saturate, and one makes a strongly cou- 
ple Color Glass Condensate. These features may re- 
semble low energy strong interactions more than high 
energy, and for the reasons stated is difficult to apply 
to jet quenching . [84] - {86] 

Nevertheless, the experimental data on jet energy 
loss and heavy particle flow is stunning. This is par- 
ticularly so for the distribution of particles produced 
by the backward going partner of a tagged jet.. The 
distribution of such particles is double peaked around 
180 degrees opposite in azimuthal angle of the tagged 
jet. This has led to speculations that the effect may 
arise from a Mach cone, that is the backward going jet 
propagating supersonically in the QGP media. Such 
Mach cones are clearly seen in the AdSCFT simula- 
tions, Fig. 



14 



9. Quarkyonic Matter 

Matter at high baryon density was the subject of 
studies at the CERN SPS and the AGS. There will 
soon be a study using low energy beams at RHIC. 
There have been a number of speculations concerning 
possible Color Superconductivity of this matter, that 
may occur at very high baryon densities and quite 
low temperature. [91]- [52] Recent arguments have sug- 
gested that the phase structure of QCD matter is even 
richer than had previously been envisioned. [93 - L95J 

QCD matter was originally thought to be divided 
into a confined word of baryons and mesons and a de- 
confined world of quarks and gluons. The new fea- 
ture on this phase diagram is Quarkyonic Matter. 
This matter is confined, but the Fermi sea of quarks 
has interactions that can be treated perturbatively. 
Baryon excitation near the Fermis surface, and all 
mesonic and glueball excitations are confined. The 
word quarkyonic is a hybrid of quark and baryonic for 
this reason. 

The reason why quarkyonic matter is confined is 
that at low temperature, the only particles that can 
short out the quark potential are quark loops. In the 
large N c limit of QCD, these loops are suppressed by 
a factor of 1/N C . So at least at large N c , there is a 
world that has densities very high compared to the 
QCD scale Aq CD that is confined. 

In so far as the large N c limit is valid, the transition 
to the deconfined world should be indepedant of den- 
sity. This is confirmed in lattice computations where 
the deconfincment temperature depends only weakly 
on baryon density. 

The phase transitions in QCD are largely charac- 
terized by changes in the number of degrees of free- 
dom. The confined world has light mass pion degrees 
of freedom. The decconfined world has 2(N 2 — l) gluon 
and 4N c Np quark-antiquark degrees of freedom. The 
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Figure 14: A Mach cone generated by a fast moving particle in an AdSCFT computation 




Figure 15: The phases of QCD at finite temperature and 
density. 



quarkyonic world has the degrees of freedom of 2N c Nf 
quarks plus light mass meson degrees of freedom. 

The chiral properties of the quarkyonic world are 
not yet understood although there have been some in- 
teresting and strongly motivated conjectures [96 -[98 . 

A prominent feature of the QCD phase diagram is 
a triple point where the quarkyonic, confined and de- 
confined world coexist. There may be a nearby critical 
point if the transitions between these different phases 
for some range of temperature and density are truly 
phase transitions, as opposed to cross overs. If there is 
no true phase transition, then one expects rapid cross 
over between the various phases. 



10. Conclusions 

The strongly interacting Quark Gluon Plasma is 
well described by perfect fluid hydrodynamic compu- 
tations at RHIC. The origin of this behaviour, and jet 



quenching within the Quark Gluon Plasma has proven 
difficult to understand. 

The Color Glass Condensate are suggested by the 
data from RHIC and Hera but their is not yet consen- 
sus that this is the origin of the various phenomena 
elaborated in the paper. I believe a compelling case 
can be made after the preliminary data on the ridge 
and long range correlations is published, and data on 
forward-backward correlations at forward rapidiities 
is represented and published. 

Mater at high baryon density is less understood. 
This is a difficult issue to treat theoretically since 
numerical Monte Carlo methods that conclusively 
demonstrated the properties of matter at high tem- 
perature do not work well at high baryon density. Ex- 
perimental data on the high baryon density region is 
difficult since this involves experiments at lower en- 
ergy, and it is not clear how far one can probe the 
high density region. 

Please note that this paper contains a review of re- 
cent theoretical developments in the theory of high 
energy density strongly interacting matter. It is not 
complete in its treatment, and in particular the discus- 
sion of more established ideas in the field. It reflects 
what I consider to be of interest for many theorists 
in the field. As such, and since this is a conference 
talk, not an extensive review, the references are not 
as comprehensive as I might have preferred. 
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